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This  paper  discusses  physical  process  of  thermal  transfer  in  energy  pile  foundations,  which  function  as 
structural  support  for  the  buildings  as  well  as  heat  exchangers  serving  energy  to  the  buildings.  Derivation 
of  conservation  of  energy  balance  is  presented,  depends  on  its  type  of  heat  transfer  within  a  whole 
system  of  soil  and  foundation.  In  order  to  understand  the  mechanical  implication  behind  their  application, 
simulation  numeric  with  finite  difference  method  is  conducted,  concerning  an  energy  pile  in  homogenous 
soil  under  static  thermal  load.  The  study  takes  into  account  two  different  conditions  of  contact  between 
soil  and  pile:  perfectly  contact  and  sliding  contact  using  frictional  interface  elements.  The  results  show 
that  temperature-induced  mechanical  behaviour  of  pile  and  soil  is  strongly  related  to  the  condition  of 
contact  between  them.  Further  work  is  projected  to  consider  a  more  appropriate  law  that  corresponds  to 
cyclic  thermal  loading  of  energy  piles  due  to  its  seasonal  cooling  and  heating  operation  throughout  the 
year. 
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1.  Introduction 

Utilisation  of  shallow  geothermal  energy  in  providing  thermal 
needs  of  building  is  already  common.  Shallow  geothermal  energy  is 
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the  average  ground  temperature  at  shallow  depth  10-50  m,  where 
its  value  ranges  from  10°C  to  15°C  in  most  European  countries 
[1  ].  This  kind  of  technology  is  well  known  as  conventional  ground 
heat  exchanger,  where  ground  thermal  energy  is  extracted  from  the 
soil  via  heat  pump  connection  to  the  building.  The  ground  operates 
as  heat  source  supplying  warm  energy  to  building  during  winter 
season,  whilst  during  summer  season  it  functions  as  heat  sink  when 
cold  energy  of  building  is  required. 

There  are  two  types  of  conventional  ground  heat  exchangers. 
The  first  one  is  open  loop  system  where  the  groundwater  is  used 
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Fig.  1.  Energy  pile  components. 


as  a  heat  carrier  and  pumped  directly  from  an  aquifer  to  the  heat 
pump.  The  remained  waste  water  is  returned  to  the  same  aquifer 
or  to  another  well.  This  system  is  simpler  but  it  induces  some  envi¬ 
ronment  problems  like  groundwater  pollution  and  needs  high  cost 
of  maintenance  due  to  clogging  in  the  wells  and  heat  exchangers 
components. 

The  second  one  is  closed  loop  system  where  closed-coils 
absorber  pipes  are  embedded  into  borehole,  laid  either  horizontally 
or  vertically.  A  heat  carrier  medium  is  circulated  through  absorber 
pipes  to  transport  heat  from  the  ground  to  the  heat  pump.  This  sys¬ 
tem  is  also  called  as  borehole  heat  exchanger  (BHE)  system.  Because 
of  its  high  cost  of  drilling  (for  vertical  loop)  and  its  need  of  huge  land 
area  (for  horizontal  loop),  closed  loop  system  has  recently  been 
merged  into  structural  foundation  elements. 

Energy  foundations  have  both  functions  as  structural  elements 
and  heat  exchanger  elements.  This  system  is  like  closed-loop  con¬ 
ventional  ground  heat  exchanger  system  unless  the  absorber  pipes 
are  installed  on  foundation  elements.  The  main  advantage  of  this 
system  is  the  foundation  elements  that  serve  as  heat  exchang¬ 
ers  are  already  required  for  structural  reason  and  no  need  to 
be  constructed  separately.  Furthermore,  concrete  as  heat  transfer 
medium  in  energy  foundations  produces  higher  thermal  energy 
due  to  its  higher  thermal  capacity  than  grout  or  rock  in  bore¬ 
hole  heat  exchanger.  Therefore,  as  well  as  reduce  significantly  cost 
installation  and  land  use,  energy  foundations  are  able  to  increase 
thermal  productivity  while  preserving  environmental  sustainabil¬ 
ity. 

Some  applications  of  energy  foundations  technology  have  been 
developed  in  last  decade,  such  as:  Main  Tower  in  Frankfurt- 
Germany  on  1 999  (112  energy  piles)  [2],  Keble  College  in  Oxford-UK 
on  2001  (90  energy  piles,  74MWh  annual  heating  and  55MWh 
annual  cooling)  [3],  Dock  Midfield  Zurich  Terminal  Airport- 
Switzerland  on  2003  (306  energy  piles,  2210  MWh  annual  heating 
and  11 70  MWh  annual  cooling)  [4],  and  Lainzer  Tunnel  Vienna- 
Austria  on  2004  (59  energy  piles,  214  MWh  annual  heating)  [1].  All 
of  those  projects  confirm  that  utilisation  of  energy  foundation  sys¬ 
tem  can  save  the  annual  overall  costs,  even  though  in  the  beginning 
the  investment  costs  are  highly  expensive.  Energy  piles  applica¬ 
tion  at  Zurich  Airport  reported  cost  saving  up  to  €36,900  relative 
to  conventional  system  with  annual  thermal  cost  €0.04/kWh  [4]. 

Despite  of  its  successful  implementation,  there  are  two  main 
problems  of  this  technology  to  be  considered.  One  concerns  about 


sustainability  of  thermal  capacity  in  the  system,  another  concerns 
about  mechanical  durability  of  structure  foundation  due  to  thermal 
solicitation  [5].  Derived  from  those  main  problems,  this  paper  will 
discuss  separately  each  problem  of  particularly  energy  piles  into: 

(a)  physical  process  of  heat  transfer  which  is  influenced  by  seasonal 
ground  temperature  and  groundwater  flows  (Section  3), 

(b)  thermal  transfer  balance  between  soil-concrete  foundations- 
and  building  (Section  4), 

(c)  mechanical  implication  induced  by  additional  thermal  loads 
which  change  the  behaviour  of  structure  (Section  5). 

In  order  to  understand  energy  piles’  thermo-mechanical 
behaviour,  preliminary  study  of  an  energy  pile  under  static  ther¬ 
mal  loads  in  each  operation  season  is  being  simulated.  Since  the 
response  of  piles  to  vertical  loads  strongly  depends  on  the  char¬ 
acteristics  at  the  soil-pile  interface,  it  is  of  major  interest  to  take 
into  consideration  the  interface  condition  in  the  analysis  of  the 
behaviour  of  piles  under  thermal  loading  [6,7].  In  this  study,  two 
comparisons  of  contact  conditions  are  proposed:  perfect  contact 
and  sliding  contact  between  soil  and  pile  using  frictional  interface 
elements  (Section  6). 

2.  Energy  piles  system 

Energy  piles  system  contains  of  three  principal  circuits  (Fig.  1). 
The  first  circuit  is  the  soil  as  heat  source  in  winter  season  where 
ground  energy  is  withdrawn  for  heating  needs  of  building.  It 
becomes  inversely  heat  sink  in  summer  season  where  heat  is 
rejected  back  to  soil  in  order  to  recharge  ground  energy  and  to  cool 
the  building. 

The  second  circuit  contains  of  closed  absorber  pipes  in  concrete 
piles  through  which  heat  carrier  fluids  are  circulated  to  exchange 
energy  from/to  soil  to/from  building.  Heat  carrier  fluids  are  water 
with  antifreeze  solution  which  is  able  to  avoid  freezing  at  incli¬ 
nation  angle  of  absorber  pipes.  Hence,  glycol-water  mixtures  are 
the  most  suitable  solution,  containing  also  additives  to  prevent 
corrosion  in  the  header  block  of  valves  of  the  heat  pump  [8].  Con¬ 
cerning  the  type  of  absorber  pipes  used  in  energy  piles,  Hamada 
et  al.  [9]  conducted  performance  test  of  three  kinds  closed  absorber 
pipes:  U-shaped,  double  U-shaped,  and  indirect  double-pipe  in  the 
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Fig.  2.  Scheme  of  energy  piles  system  operation:  (a)  energy  extraction  and  (b)  energy  storage. 


building  project  constructed  at  Sapporo,  Japan.  As  a  result,  the  U- 
shaped  pipe  type  was  employed  from  the  viewpoint  of  economic 
efficiency  and  workability  due  to  its  smallest  heat  rejection  rate  [9]. 

The  third  circuit  is  closed-pipelines  heating/cooling  network 
which  embedded  in  the  floors  and  walls  of  building.  Heat  pumps 
are  linked  between  second  and  third  circuit  in  order  to  raise  the 
original  non-usable  heat  resources  to  a  higher  usable  room  temper¬ 
ature.  Quantity  of  energy  production  depends  on  the  efficiency  of 
heat  pumps,  which  is  described  by  the  coefficient  of  performance 
(COP).  COP  is  the  ratio  of  the  output  energy  divided  by  the  input 
energy  (electricity  for  the  compressor)  and  varies  from  3.0  to  6.0 
with  present  equipment  (the  higher  the  number  the  better  the  effi¬ 
ciency)  [10].  For  a  best  energy  produced  with  economical  cost,  heat 
pumps  require  a  number  of  COP  higher  than  4.0  [8]. 

Two  operation  schemes  of  energy  piles  are  possible:  single  mode 
of  energy  extraction/energy  input  and  seasonal  mode  operation 
with  heating  and  cooling  storage  (Fig.  2).  For  single  mode  opera¬ 
tion,  the  energy  flow  takes  place  in  one  direction  only,  for  example: 
heating  purposes  in  glacial  countries  or  cooling  purposes  in  tropical 
countries.  If  only  heating  or  cooling  is  performed,  high  perme¬ 
ability  ground  and  high  hydraulic  gradient  of  groundwater  are  of 
advantage  [8].  However,  the  seasonal  operation  uses  the  thermo¬ 
dynamic  inertia  of  the  soil  in  order  to  store  thermal  energy  in 
the  ground  for  later  operation  with  reversed  energy  flow.  Conse¬ 
quently,  the  seasonal  operation  can  produce  energy  equilibrium 
automatically  in  the  ground  over  a  complete  heating/cooling  period 
of  a  year  [1].  This  operation  system  is  preferred  by  the  ecologist 
due  to  its  environmental  friendly  for  the  sustainability  of  ground- 
water. 


Heat  transfer  in  energy  piles  system  is  surrounded  by  complex 
process.  In  the  first  circuit,  the  transfer  depends  on  soil  ther¬ 
mal  and  hydrology  properties,  which  are  specifically  conducted  by 
thermal  conductivity,  specific  heat  extraction,  and  water  content 
of  soil.  In  the  second  circuit,  thermal  transfer  affects  mechanical 
behaviour  at  concrete  piles,  predominantly  related  to  the  inter¬ 
action  between  soil-concrete  piles-and  upper  structure.  Finally, 
coefficient  of  performance  of  heat  pumps  is  the  key  to  avoid 
the  energy  loss  between  the  second  and  third  circuit.  Because  of 
the  fact  that  no  design  method  is  yet  available  to  consider  the 
complex  interaction  between  thermal  storage  and  the  mechanical 
behaviour  of  energy  piles  [11],  the  present  work  offers  some  steps 
to  implement  energy  piles’  technology  as  explained  at  flowchart  in 
Fig.  3. 

3.  Physical  process  in  heat  transfer 

Soil  is  a  porous  medium  formed  by  solid  grains,  water,  and 
gases.  Due  to  their  compositions,  heat  transfer  in  soil  consist  of 
many  mechanism,  including:  conduction  through  the  soil  grains, 
liquid,  and  gases;  convection  through  liquid  and  vapour  diffusion; 
latent  heat  transfer  through  evapouration-condensation  cycles 
during  the  phase  change  of  water;  radiation  in  gas-filled  pores; 
and  freezing-thawing  process  that  should  be  avoided  in  energy 
foundations  system  [8,12].  Heat  conduction  leads  dominantly,  heat 
convection  in  liquid  diffusion  occurs  by  presence  of  groundwater 
flow,  vapour  diffusion  and  latent  heat  transfer  might  only  possi¬ 
ble  in  unsaturated  soil,  and  radiation  contributes  negligibly  to  heat 
transfer. 
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Fig.  3.  Flow  chart  of  design  steps. 


3.1.  Ground  temperature 

Burger  et  al.  [  1 3  ]  divides  the  soil  into  three  zones  of  temperature : 
heterothermal,  neutral,  and  homothermal  zone.  External  temper¬ 
ature  and  radiation  affect  ground  temperature  in  heterothermal 
zone.  In  this  zone,  temperature  varies  in  time  and  depth  as  function 
T(z,t):  where  T(t)  defines  thermal  cycle  loads  in  annual  exter¬ 
nal  temperature  and  T(z)  indicates  thermal  waves’  propagation  in 
depth.  This  is  the  most  critical  zone  depth  that  should  be  properly 
taken  in  account  of  energy  piles  design  calculation.  The  function  of 
temperature  forms  a  sinusoidal  wave,  as  following  equation: 

T(z,  t)  =  rave  +  A0e~{z/d)  f sin  f wt  -  ^  ^  ( 1 ) 

where  raVe  is  the  average  annual  soil  temperature;  A0  is  maximum 
annual  amplitude;  z  is  the  depth  in  which  temperature  is  investi¬ 
gated;  d  is  the  damping  depth  of  annual  fluctuation;  and  co  is  the 
annual  radial  frequency. 

Fluctuation  of  temperature  becomes  constant  since  it  arrives 
in  neutral  zone  at  depth  10m  until  50m  below  the  surface  [13] 
and  remains  closely  at  average  air  temperature,  where  the  value 
depends  on  altitude  h : 

Tair  =  11.3  -0.52/1  (2) 

In  homothermal  zone,  temperature  is  no  longer  influenced  by 
external  thermal  flux  and  no  longer  varies  by  annual  sinusoidal 


waves.  This  zone  is  where  the  deep  geothermal  or  internal  terres¬ 
trial  flux  takes  place.  Temperature  increases  regularly  in  depth  with 
geothermal  gradient  3  °C  per  100  m  [14]. 

3.2.  Soil  thermal  properties 

Heat  flux  propagation  in  soil  is  transmitted  by  heat  capacity  and 
thermal  conductivity.  Volumetric  heat  capacity  Cv  is  a  product  of 
mass  volume  p  with  specific  heat  extraction  c  that  indicates  the 
ability  of  a  substance  to  store  heat  energy.  The  greater  its  heat 
capacity,  the  more  heat  it  can  gain  or  lose  per  unit  rise  or  fall  in 
temperature.  Whereas  thermal  conductivity  X  is  the  quantity  of 
heat  transferred  through  a  unit  area  of  the  conducting  body  in  unit 
time  under  a  unit  temperature  gradient  [15].  This  indicates  the  rate 
at  which  heat  will  be  transferred  between  exchanger  loops  and 
the  surrounding  transfer  medium  for  a  given  temperature  gradi¬ 
ent.  Thermal  conductivity  of  the  surrounding  transfer  medium  is 
the  critical  value  that  determines  the  number  of  heat  exchanger 
pipes  required  to  produce  energy  needs  of  building  by  pumping 
specific  fluid  temperature  through  the  ground  loops. 

Thermal  properties  of  soil,  including  heat  capacity  and  thermal 
conductivity,  cannot  be  separated  with  mineral  composition  of  soil 
itself  as  mentioned  in  Eqs.  (3)  and  (4)  with  index  “s”  denotes  solid 
composition  and  denotes  fluid  composition.  Porosity  factor  n, 
water  and  air  content  will  lead  the  quantity  of  heat  produced.  Moist 
or  saturated  soils  have  greater  thermal  capacity  than  dry  soils,  since 
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thermal  conductivity  for  dry  soil  is  about  five  times  lower  than  that 
of  a  saturated  soil  [14]. 

Cy  —  ( 1  n)pscs  +  npfcf  (3) 

A.  =  (l  -n)Xs  +  nXf  (4) 

3.3.  Groundwater  flow 


As  water  and  air  content  take  important  place  in  producing  ther¬ 
mal  capacity  of  soil,  the  presence  of  groundwater  flow  is  one  of  key 
factor  in  heat  transfer  mechanism.  In  fully  saturated  soil  with  no 
presence  of  pore  gas,  groundwater  composes  only  of  liquid  flow, 
so  that  heat  convection  occurs  in  liquid  diffusion.  In  unsaturated 
soil  where  the  pore  gas  presences  inevitable,  groundwater  flow  has 
both  of  liquid  and  vapour  flow,  so  that  heat  convection  counts  both 
liquid  and  vapour  diffusion.  It  is  also  possible  to  have  a  water  phase 
change  into  fully  vapour  that  correlates  with  vapour  latent  transfer. 
A  widely  used  theory  for  coupled  heat  and  moisture  flow  through 
soils  was  developed  by  Philip  and  de  Vries  [16]  for  both  liquid  and 
vapour  phase  flows. 

Groundwater  flow  in  soil  varies  under  water  content  0  and  tem¬ 
perature  T  gradients.  Water  content  variation  is  led  by  isothermal 
diffusivity  De  while  temperature  variation  is  controlled  by  thermal 
diffusivity  DT.  In  general,  water  content  gradient  is  less  than  those 
of  temperature,  therefore  the  flux  due  to  water  content  gradient 
is  more  important  than  those  due  to  temperature  gradient,  even 
with  the  maximum  temperature  gradient  [16].  Hadley  and  Eisen- 
stadt  [17]  have  observed  that  moisture  transfer  under  temperature 
gradient  is  negligibly  small  both  in  very  dry  (unsaturated)  and  in 
very  wet  (saturated)  media  [16,17]. 

In  liquid  phase,  the  velocity  of  groundwater  is  given  by  Eqs. 
(5)-(7),  referring  to  Darcy’s  equation  where  k  is  hydraulic  con¬ 
ductivity  and  i  is  unit  vector  in  vertical  direction.  Here,  isothermal 
liquid  diffusivity  Dqj  and  thermal  liquid  diffusivity  Dr/  depend  on 
capillary  suction  head  xj/  and  surface  capillary  tension  a. 


V[  =  -De jgrad  9  -  DT  igrad  T  -  ki 

(5) 

0J  ~  k  30 

(6) 

Dtj  ~k0dT 

(7) 

Movement  in  vapour  phase  is  a  process  of  the  diffusion  of  water 
vapour  in  the  air-filled  pores  which  can  be  approximated  by  mod¬ 
ifying  the  Fields  law  of  diffusion: 


vv  =  -Dj  vgrad  T  -  Dd  vgrad  0 


where  DT  v  is  the  thermal  vapour  diffusivity  and  De  v  is  the  isother¬ 
mal  vapour  diffusivity. 

Groundwater  flow  is  then  a  sum  of  flow  in  liquid  and  vapour 
phase.  The  two  liquid  diffusivities  tend  to  be  the  most  important 
ones  at  high  moisture  contents,  whilst  the  two  vapour  diffusivities 
are  dominant  at  low  moisture  contents  [16].  Refers  to  precedent 
hypotheses,  groundwater  flow  in  saturated  soil  is  given  in  Eq.  (9) 
while  in  unsaturated  soil  is  mentioned  in  Eq.  (10). 


Vsat  =  -kytgrad  0  = -k  grad  f 


Vunsat  —  -De  vgrad  0 


(10) 


Site  installation  of  Finkernweg  at  Switzerland  has  been  con¬ 
ducted  by  Fromentin  et  al.  [5],  concerning  75  energy  piles  with  11m 
active  length  to  produce  90  kW  and  7000  m3  heat  energy.  Thermal 
measurement  is  investigated  during  February  until  May  1995.  It  is 
shown  that  groundwater  flow  does  not  have  significant  influence 
in  short  term  operation.  However,  groundwater  flow  may  repro¬ 
duce  heat  concentration  energy  around  the  piles  with  velocity  up 


to  lOcm/day  in  long-term  period  [5]  so  that  it  plays  an  important 
role  for  pile’s  dimensioning.  Also,  if  a  flow  of  water  is  present  and 
sufficiently  large,  a  natural  regeneration  of  the  soil  is  achieved  [14]. 
As  a  result,  the  amount  of  heat  extracted  during  winter  (for  heating) 
no  longer  depends  on  the  amount  injected  in  the  summer  (for  the 
cooling).  Hence,  groundwater  flow  is  of  advantage  to  single  mode 
energy  extraction  system  but  further  study  is  required  to  ensure 
sustainability  of  ground  water. 

3.4.  Heat  flux 

The  rate  of  energy  transfer  in  soil  through  solid,  liquid,  or  gas 
pores  is  called  as  heat  flux.  It  is  important  to  define  the  heat  flux  in 
order  to  quantify  the  energy  produced. 

For  heat  conduction  mechanism,  heat  flux  is  based  on  Fourier 
law: 

Qcond  =  ~k  grad  T  (11) 

while  in  liquid  convection  mechanism  is  determined  by: 


Qconv  —  PfCfVl(T  —  To)  (12) 

It  should  be  noted  that  groundwater  flow  in  heat  convection 
flux  is  the  groundwater  flow  in  liquid  diffusion.  In  reverse  when 
we  consider  latent  heat  flux,  groundwater  flow  which  contributes 
in  energy  transfer  is  groundwater  flow  in  vapour  diffusion,  as  men¬ 
tioned  as: 

Qlatent  =  LyVy  (13) 


3.5.  Energy  conservation 


Finally,  heat  transfer  in  soil  should  satisfy  the  energy  conserva¬ 
tion  of  thermodynamics  law,  which  is: 


c'¥+v’"° 


(14) 


Table  1  summarises  heat  transfer  mechanism  for  each  type  of 
soils  and  their  energy  conservations  are  given  in  Eqs.  ( 1 5)— ( 18). 
For  saturated  soil  with  no  groundwater  flow: 


Cv—=div(\  grad  T), 
ot 

for  saturated  soil  with  groundwater  flow: 
dT 

Cv— -  pfCfk  grad  xj/  grad  T  =  div(X  grad  T), 

for  unsaturated  soil  without  groundwater  flow: 
dT 

Cv—=div(X  grad  T)-Lv  div{D0vap  grad  0), 

and  for  unsaturated  soil  with  groundwater  flow: 
dT 

Cv— -  pfCfk  grad  xfr  grad  T  =  div(X  grad  T) 


Ly  div{De  vap  grad  0) 


(15) 


(16) 


(17) 


(18) 


4.  Thermal  transfer  balance 

Energy  piles  are  dual-purpose  as  structural  foundation  and  heat 
exchanger.  Thermal  capacity  of  their  heat  exchanger  elements 
will  take  a  role  for  structural  dimensioning  of  their  foundation 
elements.  Their  functions  to  satisfy  thermal  and  mechanical  (struc¬ 
tural)  capacity  should  be  carefully  designed  by  mechanical  and 
geotechnical  engineers.  For  Q3  is  the  energy  required  of  building  at 
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Table  1 

Division  of  heat  transfer  process  in  soils. 
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Heat  transfer 

Saturated 

Unsaturated 

No  ground-water  flow 

With  ground-water  flow 

No  ground-water  flow 

With  ground-water  flow 

In  solid  conduction 

V 

V 

V 

V 

In  liquid  convection 

X 

V 

X 

V 

In  vapour  latent 

X 

X 

V 

V 

Table  2 

Properties  of  model. 

Thermo-elastic  properties 

Unit 

Soil 

Concrete 

Density 

P 

Pa 

1950 

2500 

Bulk  modulus 

E 

MPa 

10 

20,000 

Shear  modulus 

G 

MPa 

3.75 

7500 

Thermal  conductivity 

X 

W/m2 

1.5 

1.8 

Specific  heat  extraction 

c 

J/kg°C 

800 

880 

Coefficient  of  thermal  expansion 

a 

J/°c 

5  x  10-6 

1  x  10“5 

Thermo-elastoplastic  properties 

Unit 

Interface 

Normal  stiffness  kn 

MN/m 

22 

Shear  stiffness  ks 

MN/m 

8.33 

Cohesion  c 

kPa 

1 

Friction  angle  <fi 

0 

30 

third  circuit  and  Q2  is  input  energy  at  heat  pump  at  second  circuit, 
the  energy  conservation  which  occurs  in  the  heat  pump  is: 


Q3  =  Q2  +  P HP 

M'-asO-42 


(19) 

(20) 


Input  energy  at  heat  pump  Q2  is  equal  to  resources  heat  energy 
from  soil  through  concrete  piles  Qi.  Hence,  to  design  the  number 
n  of  concrete  piles  required,  Eqs.  (21)  and  (22)  are  given  as  below, 
where  l  is  the  length  of  material. 

02  =  Qi  (21) 


J  ( nqCvdl ) 


concrete  ~ 


( QCpdl^soU 


(22) 


To  minimise  energy  loss  between  each  circuit,  high-raised 
energy  can  be  produced  by  the  help  of  heat  pumps  performance. 
Study  of  thermal  performance  of  energy  piles  applied  at  Zurich  Air¬ 
port  showed  that  annual  heat  pump  performance  coefficient  is  3.9. 
The  overall  system  efficiency,  defined  by  the  ratio  of  the  thermal 
energy  delivered  by  the  system  (heating  and  cooling)  over  the  total 
electric  energy  required  to  run  it,  is  established  to  5.1  [4].  With 
a  high  efficiency  system,  energy  piles  will  maintain  in  long-term 
operation  with  less  energy  loss. 


5.  Mechanical  equilibrium 

In  general,  the  mechanical  response  of  a  system  must  sat¬ 
isfy  equilibrium,  compatibility,  material  constitutive  behavior  and 
boundaries  conditions  [18].  Among  these  conditions,  thermal  effect 
is  involved  in  mechanical  response  of  energy  piles.  Indeed,  the 
change  of  temperature  both  in  concrete  and  soil  will  produce 
supplementary  thermal  deformation  that  will  modify  the  strain 
state  of  materials.  The  quantity  of  thermal  deformation  £th  is  a 
product  of  thermal  expansion  coefficient  a  with  gradient  temper¬ 
ature  occurred  as  mentioned  in  Eq.  (23).  Elastic  deformation  se 
will  change  in  function  of  thermal  deformation  obtained  and  total 
deformation  s  caused  by  static  weights  of  system  structure,  given 
in  Eq.  (24). 


(23) 


£e  =  £  -  £th  (24) 

Since  lateral  movements  of  concrete  pile  are  restrained  by  soils, 
thermal  deformation  will  mainly  occurs  in  vertical  axis  of  pile.  The 
vertical  displacements  are  not  uniform  in  depth  and  strongly  relates 
to  lateral  friction  mobilised  in  the  shaft  of  pile,  in  contact  with  soil 
[19].  The  supplementary  thermal  deformation  brings  out  increas¬ 
ing  of  mobilised  shaft  friction  at  interface  soil-concrete  piles  also 
uplift  displacement  and  additional  thermal  compressive  stress  at 
interface  concrete  piles-upper  structure  [20].  The  relation  between 
axial  force  Q,  mobilised  shaft  friction  qs ,  and  deformation  at  a  given 
depth  z  is  expressed  on  following  equation  [21  ]: 


Q(z)  =  £e{z)AEpile 


(25) 


Qs(z) 


Q.(z) 

Mz) 


AEpiie  A £e(z) 
As(z)  dz 


(26) 


where  Epue  is  young  modulus  of  concrete  pile,  A  is  area  of  pile’s 
section,  As  is  shaft  area  at  certain  depth. 

A  four  storey  building  built  with  97  piles  was  carried  out 
to  an  in  situ  test  of  a  heat  exchanger  pile  at  the  Swiss  Federal 
Institute  of  Technology  in  Lausanne  (EPFL).  Data  from  tests  indi¬ 
cate  that  a  temperature  increment  of  1  °C  results  in  an  additional 
temperature-induced  vertical  force  on  the  order  of  100 1<N  [22].  As 
a  consequence,  the  total  axial  load  in  the  pile  is  twice  as  large  as  the 
one  due  to  purely  mechanical  loading,  with  a  large  solicitation  of 
the  toe  because  the  pile  toe  sits  on  very  stiff  sandstone  that  restrain 
the  pile. 

Another  in  situ  test  of  energy  piles  was  conducted  in  Lambeth 
College,  United  Kingdom  with  a  total  143  piles.  The  different  with 
those  from  Lausanne  test  is  in  this  test,  the  piles  are  floating  in  the 
ground  so  that  both  ends  of  pile  are  free  to  move.  When  a  cooling 
cycle  is  applied  to  the  pile,  the  pile  contracts  and  leads  to  induce 
tensile  force.  Along  the  upper  part  of  the  shaft,  mobilised  lateral 
friction  on  the  interface  soil-pile  will  increase  in  the  same  sense 
as  that  mobilised  by  compression  loading  applied  at  the  pile  head, 
while  in  the  lower  part  of  the  pile  will  be  in  the  opposite  sense 
so  reduce  the  mobilised  friction  [23].  The  reverse  responses  are 
appeared  when  a  heating  cycle  is  applied,  which  suggest  the  pile 
behaves  in  purely  thermo-elastic. 

If  the  energy  piles  work  on  seasonal  mode  operation  sys¬ 
tem,  they  are  subjected  by  alternate  thermal  loads.  The  alternate 


£th  =  a  AT 
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Fig.  4.  3D  mesh  discretisation  of  model:  (a)  one-fourth  of  model  and  (b)  interfaces  location  at  Model  2. 


contraction/expansion  piles  create  cyclic  behaviour  at  interface 
soil-concrete  piles-and  upper  structure.  As  explained  in  Section 
2,  from  an  energy  balance  point  of  view,  seasonal  mode  operation 
system  is  the  most  favourable  one  referring  to  energy  balance  point 
of  view.  But  if  we  look  from  the  view  of  mechanical  durability, 
cyclic  behaviour  at  interface  will  decrease  strength  resistance  of 
structure. 

6.  Case  study:  temperature-induced  mechanical  behaviour 
of  an  energy  pile  under  static  thermal  loading 

6.2.  Numerical  model 

Numerical  method  such  as  finite  difference  method  provides  an 
efficient  tool  for  solving  complicated  partial  differential  equations 
such  as  those  that  govern  the  thermo-hydro-mechanical  problem 
of  energy  pile,  because  of  its  ability  to  integrate  complexities  aris¬ 
ing  from:  geometry,  thermo-hydro-mechanic  loading,  interaction 
between  different  rigidity  bodies,  etc.  The  case  under  study  con¬ 
cerns  with  the  analysis  of  an  energy  pile  under  static  thermal 
loading.  A  single  pile  with  a  square  section  of  B  =  60  cm  and  length 
L=  15  m  is  founded  into  a  homogeneous  soil  mass.  Both  of  pile  and 
soil  are  assumed  to  behave  into  linear  thermo-elastic  domain.  Prop¬ 
erties  of  materials  are  shown  in  Table  2,  where  thermal  properties 
are  chosen  in  reference  to  Hillel  [15]. 

Analyses  are  carried  out  using  finite  difference  code  FLAC30 
[24].  Fig.  4a  shows  the  numerical  mesh  used  into  the  numerical 
modelling.  Owing  to  the  symmetry  of  the  problem  (symmetry  of 
load  and  geometry),  only  one-fourth  of  the  complete  domain  is 
modelled.  The  mesh  was  fixed  after  a  parametric  analysis  con¬ 
cerning  the  lateral  extension  of  the  soil  domain.  The  soil  is  refined 
around  the  pile,  in  order  to  increase  precision  in  the  areas  of  high 
strain  gradient.  According  to  Karthigeyan  et  al.  [25]  and  Mroueh 
and  Shahrour  [7],  the  soil  mass  dimension  depends  on  the  pile 


section  and  length.  Lateral  extension  of  the  domain  is  equal  to  1 5  m 
(25B).  The  height  of  soil  mass  is  equal  to  30  m  (2  L).  As  mentioned 
in  Section  1 ,  the  response  of  piles  to  vertical  loads  strongly  depends 
on  the  characteristics  at  the  soil-pile  interface  [6,7].  This  study 

T(°C) 


0  5  10  15  20  25  30  35 


- Spring  - Summer - Autumn  - Winter 

Fig.  5.  Fluctuation  of  initial  ground  temperature. 
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Fig.  6.  Initial  condition  of  ground  temperature. 


proposes  comparison  of  two  conditions  of  soil-pile  contact  in 
order  to  understand  their  different  behaviour.  Model  1  considers 
a  perfect  contact  between  soil  and  energy  pile  that  means  there  is 
no  relative  movement  between  the  nodes  of  pile  and  soil.  Model  2 
introduces  frictional  model  for  the  soil-pile  interface,  governed  by 
an  elastic  perfectly  plastic  behaviour.  These  interfaces  elements 
are  able  to  localise  relative  movements  and  frictions  mobilised 
between  two  different  solids  (Fig.  4b).  Elasticity  part  is  described 
by  the  following  relation: 


where  r  and  an  denote  the  shear  and  normal  stress  and  st  and  sn 
indicate  the  shear  and  normal  strain  at  the  interface,  respectively. 
The  interface  shear  stiffness  ks  and  the  normal  stiffness  kn  are 
chosen  in  accordance  with  interfaces’  theory  and  background  of 
the  code  FLAC30  [24].  Plasticity  part  is  defined  by  a  yield  surface 
based  on  the  Coulomb  failure  criterion: 

\r\  =  crn  tan  (pi  +  C[  (28) 

where  (pt  and  q  are  the  friction  angle  of  interface  and  cohesion 
factor  at  interface.  The  interfaces  properties  are  shown  in  Table  2. 

At  first,  soil  system  is  initialised  by  sinusoidal  surface  tempera¬ 
ture  with  average  ground  temperature  14  °C.  The  soil  has  internal 
terrestrial  heat  volume  r  as  0.001  W/m3  and  terrestrial  surface  flux 


ks  0 
0  kn 


(27) 


Q(kN) 

0  5  10  15  20  25  30 


Fig.  7.  Normal  forces  due  to  gravity  load. 


q  at  depth  30  m  as  0.0544  W/m2 .  Due  to  symmetrical  condition,  flux 
null  is  imposed  in  symmetrical  axes.  The  presence  of  groundwater 
flow  is  neglected  in  assuming  that  soil  is  fully  saturated  so  that  heat 
diffusion  works  in  purely  conduction.  This  sinusoidal  surface  tem¬ 
perature  produces  fluctuation  of  ground  temperature  in  depth  and 
in  time  that  becomes  constant  near  to  average  ground  temperature 
at  the  depth  >1 0  m  (Fig.  5).  This  is  justified  by  theoretical  hypothesis 
in  Section  3.1. 

According  to  that  fluctuation,  seasonal  inlet  fluid  temperature  is 
applied  into  the  pile  and  assumed  to  be  uniform  along  the  pile’s  axis. 
The  model  is  conducted  in  two  separated  operation  time:  (a)  winter 
mode  operation  with  inlet  temperature  is  set  on  5  °C  (cooling  cycle 
-  Model  W)  and  (b)  summer  mode  operation  with  inlet  temperature 
25  °C  (heating  cycle  -  Model  S).  Initial  temperature  of  soil  before 
the  injection  of  inlet  fluid  at  each  season  is  shown  in  Fig.  6. 

The  pile  is  only  subjected  by  gravity  load  and  then  by  temper¬ 
ature  applied  during  six  months  loading  for  each  operation  mode. 
No  service  load  coming  from  upper  structure  is  considered  to  prop¬ 
erly  understand  the  effect  of  thermal  load  in  pile’s  mechanical 
behaviour.  Fig.  7  shows  the  initial  normal  forces  due  to  gravity  load 
in  two  different  conditions  of  contact.  For  the  clarity,  positive  sign 
signifies  compression  on  forces  and  settlement  on  displacements. 
Model  2  with  interfaces  imposed  in  the  zone  of  contact  gives  reduc¬ 
tion  of  normal  forces  in  pile  because  interface  zones  take  a  role 
to  transfer  the  forces  more  uniformly  to  the  system  by  reducing 
discontinuity  of  residual  forces. 

Comparison  of  thermal  diffusion  by  the  function  of  contact’s 
condition  shows  a  remarkable  influence  (Figs.  8  and  9).  Within  the 
presence  of  interface  elements,  there  is  high  concentration  of  ther¬ 
mal  field  at  surface  contact  between  pile  and  soil.  This  leads  to  a 
statement  that  interface  behaviour  depends  also  on  temperature  as 
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Contour  of  Temperature 
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Fig.  8.  Field  of  temperature  after  cooling  cycle:  (a)  Model  1  perfect  contact  and  (b)  Model  2  with  interfaces. 
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Fig.  9.  Field  of  temperature  after  heating  cycle:  (a)  Model  1  perfect  contact  and  (b)  Model  2  with  interfaces. 


a  function  parameter,  that  interface  nodes  take  a  role  on  tempera¬ 
ture  transfer.  It  is  suggested  to  conduct  further  study  to  look  about 
interface  behaviour  law  with  the  coupling  of  temperature  function. 

Concerning  to  the  mechanical  behaviour  due  to  thermal  loading, 
thermal  displacement  of  the  surface  and  along  the  pile  are  shown 
in  Figs.  10  and  11.  Winter  mode  operation  produces  a  settlement 
at  surface  plus  shortening  of  pile  while  summer  mode  operation 


horizontal  distance  by  radius  of  pile  (R) 
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Fig.  10.  Thermal  surface  displacements. 


leads  to  an  uplift  surface  with  expansion  of  pile.  Contraction  in 
winter  mode  induces  tensile  forces  at  the  pile  while  dilatation  in 
summer  mode  adds  compressive  forces  at  the  pile.  Model  2  with 
interfaces  imposed  gives  clearly  different  mechanical  responses 
and  also  reduces  the  value  of  temperature-induced  normal  forces 
produced  (Fig.  12). 

Fig.  13  shows  frictions  mobilised  along  interface  soil-pile  at 
Model  2.  Due  to  its  contraction  in  winter  mode,  frictions  mobilised 
along  the  upper  part  of  pile  increase  whereas  along  the  lower  part 
of  pile  decrease.  Inversely  in  summer  mode,  a  reduction  of  frictions 
mobilised  is  appeared  at  upper  part  of  pile  while  the  lower  part  of 
pile  has  additional  frictions  mobilised.  At  middle  part  of  pile,  the 
friction  mobilised  tends  to  be  constant  as  what  it’s  called  null  point. 
This  behaviour  is  only  occurred  when  both  ends  of  pile  are  freely 
restrained.  These  results  are  in  good  agreement  with  in  situ  tests 
conducted  at  Lambeth  College,  UK  [23]. 


6.2.  Perspectives 

According  to  this  preliminary  numerical  model,  static  thermal 
load  in  only  cooling/heating  operation  give  slight  influence  on 
mechanical  behaviour  of  pile.  Alternate  contraction  and  dilatation 
during  seasonal  operation  time  is  predicted  to  be  a  dominant  factor 


pile's  length  pile's  length 
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wth  (mm) 


qsth(kPa) 


-1.5  -1  -0.5  0  0.5  1  1.5 


Fig.  11.  Thermal  vertical  displacements  of  pile. 
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Fig.  13.  Temperature-induced  frictions  mobilised  at  interface  soil-pile. 
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Fig.  12.  Temperature-induced  normal  forces  in  the  pile. 


that  provokes  a  strong  degradation  of  structure’s  resistance,  espe¬ 
cially  at  interface  soil-structure. 

Our  perspective  is  to  determine  the  interface  behaviour  between 
soil-structure  due  to  thermal  cyclic  solicitation.  To  achieve  this 
goal,  it  is  better  to  consider  an  appropriate  law  that  can  define  cyclic 
behaviour  of  interface  soil-structure  in  taking  into  account  plastic¬ 
ity  condition.  It  should  consider  also  the  effect  of  delta  temperature 
occurred  at  interface  behaviour.  On-going  works  have  developed 
a  constitutive  law  for  interface  soil-structure  under  mechanical 
cyclic  loading,  based  on  previous  work  of  Shahrour  and  Rezaie  [26]. 
This  law  takes  into  account  kinematic  hardening  rules  that  allow 
to  describe  kinematic  surface  in  each  cycle  loading.  The  further 
research  aims  to  adopt  this  constitutive  law  into  thermal  regime, 
coupled  with  mechanical  service  loads. 

7.  Conclusions 

Energy  piles  are  bi-function  foundations  which  support  the 
loads  of  structure  and  provide  thermal  energy  to  the  building 
as  heat  exchangers.  They’re  using  shallow  geothermal  energy  as 
heat  source,  successfully  considered  as  an  innovative  environmen¬ 
tally  friendly  building’s  technology.  Energy  piles  can  work  in  two 
system  operations:  (a)  only  cooling/heating  operation  system  for 
glacial/tropical  countries  and  (b)  seasonal  cooling-heating  opera¬ 
tion  system  for  four-season  countries.  The  latter  one  is  the  most 
favourable  for  its  automatic  energy  balance  throughout  the  year. 

This  paper  covered  a  sum  of  literature  study  to  identify  the 
remaining  problems  in  energy  piles’  application.  Quantity  of  heat 
energy  produced  depends  on  type  of  soil,  saturation  degree,  pres¬ 
ence  of  groundwater  flow,  and  of  course  heat  pump’s  specification. 
On  the  other  hand,  delta  temperature  occurred  delivers  additional 
thermal  deformation  and  thermal  stress  that  change  the  behaviour 
of  material,  especially  friction  mobilised  at  interface  soil-concrete 
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piles  also  head  uplift  displacement  and  axial  forces  at  interface 
concrete  piles-upper  structure. 

A  preliminary  numerical  simulation  is  conducted  with  finite 
difference  method  under  static  thermal  load  for  only  cooling  and 
only  heating  operation.  The  study  considered  two  different  condi¬ 
tions  of  contact  between  soil  and  an  energy  pile:  perfectly  contact 
and  sliding  contact  with  frictional  interface  elements.  In  this  case, 
the  pile  is  freely  restrained  both  at  its  ends.  The  results  show  that 
by  the  presence  of  interfaces  at  zone  of  contact  between  soil  and 
pile,  the  stresses  and  displacements  obtained  are  lower  than  those 
with  perfect  contact.  Cooling  cycle  in  winter  mode  results  in  pile’s 
contraction  so  that  produces  surface’s  settlement  and  temperature- 
induced  tensile  forces  at  pile.  If  interfaces  are  imposed  at  surface  of 
contact,  cooling  cycle  increases  frictions  mobilised  at  the  upper  part 
of  pile  while  reduces  at  the  lower  part  of  pile.  Reverses  responses 
are  obtained  at  heating  cycle  in  summer  mode.  A  higher  gradient 
temperature  applied  will  possibly  create  negative  friction  at  inter¬ 
face  soil-structure.  There  is  a  node  at  the  middle  part  of  pile  that  is 
not  influenced  by  temperature-induced  friction  mobilised,  which  is 
called  null  point.  This  behaviour  of  energy  pile  under  thermal  load 
is  none  of  difference  with  the  results  of  in  situ  tests  at  Lausanne 
and  Lambeth  College. 

Numerical  simulation  under  thermal  static  loads  for  only  cool¬ 
ing/heating  operation  pointed  out  that  there  is  slight  influence  in 
mechanical  behaviour  of  pile.  The  problem  in  this  kind  of  operation 
system  lies  in  quantity  of  energy  balance.  Mechanical  behaviour 
will  be  affected  in  seasonal  operation  time  since  the  piles  are  sub¬ 
jected  by  cyclic  thermal  loadings  during  annual  operation.  It  is 
suggested  to  concentrate  on  behaviour  of  interface  due  to  cyclic 
solicitation  in  thermal  regime  with  appropriate  cyclic  behaviour 
law  of  interface  soil-structure  considering  coupling  with  temper¬ 
ature  function  for  further  studies. 
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